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A characteristic feature of loop quantization of the isotropic and Bianchi-I space- 
times is the existence of universal bounds on the energy density and the expansion 
and shear scalars, independent of the matter content. We investigate the proper- 
ties of these physical quantities in Bianchi-II and Bianchi-IX spacetimes, which have 
been recently loop quantized using the connection operator approach. Using the ef- 
fective Hamiltonian approach, we show that for Bianchi-II spacetime, energy density 
and the expansion and shear scalars turn out to be bounded, albeit not by universal 
values. In Bianchi-IX spacetime, when the approach to the classical singularity is 
isotropic, above physical quantities are bounded. In addition, for all other cases, 
where the approach to singularities is not isotropic and effective dynamics can be 
trusted, these quantities turn out to be finite. These results stand in sharp distinc- 
tion to general relativity, where above physical quantities are generically unbounded, 
leading to the break down of geodesic equations. In contrast to the isotropic and 
Bianchi-I models, we find the role of energy conditions for Bianchi-II model and the 
inverse triad modifications for Bianchi-IX to be significant to obtain above bounds. 
These results bring out subtle physical distinctions between the quantization us- 
ing holonomies over closed loops performed for isotropic and Bianchi-I models, and 
the connection operator approach. We find that qualitative differences in physics 
exist for these quantization methods even for the isotropic models in the presence 
of spatial curvature. We highlight these important differences in the behavior of 
the expansion scalar in the holonomy based quantization and connection operator 
approach for isotropic spatially closed and open models. 



I. INTRODUCTION 

Singularities in general relativity are primarily characterized by the divergences in the 
curvature invariants and the break down of geodesic evolution. A central piece in the 
singularity theorems of Hawking, Penrose and Geroch is the Raychaudhuri equation, which 
determines the evolution of the congruence of geodesies via the properties of the expansion 
[9) and shear (cx^) scalars, and the components of the stress energy tensor. As the singularity 
is approached, these quantities blow up, resulting in the inextendability of geodesies. The 
scalars 6 and a^, also capture the extrinsic and the Weyl curvature of the spacetime, and 
hence turn out to be useful measures to gain insights on the behavior of curvature invariants 
and the nature of singularities. It is generally believed that a quantum theory of gravity 
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will shed important insights on the resolution of singularities. If such a theory allows an 
effective spacetime description to understand the behavior of geodesies, it is pertinent to 
ask whether 9 and cr^, along with the components of the stress-energy tensor are bounded 
by the quantum gravitational effects, and if yes, under what conditions. Understanding of 
these properties is vital to gain insights on the generic resolution of singularities in quantum 
gravity and the underlying requirements for geodesic completeness. 

In recent years, a lot of progress has been made in the quantization of homogeneous 
spacetimes in loop quantum cosmology (LQC) to pose such questions. LQC is a non- 
perturbative canonical quantization of homogeneous cosmologies, based on loop quantum 
gravity (LQG), which predicts resolution of singularities in various situations These 
include spatially fiat {k = 0) isotropic model sourced with a massless scalar field [2|-0], in 
presence of cosmological constant ^5ri7|] and inflationary potential [8|], spatially closed (fc = 1) 
model 10], spatially open [k = —1) model [ll|, 12 , Bianchi-I model il3i] and Bianchi-II 



Ij] and Bianchi-IX spacetimes with a massless scalar field [15|. In all of these models, the 
quantum Hamiltonian constraint turns out to be non-singular, which is a direct consequence 
of the underlying quantum discreteness predicted by LQG. Further, one recovers classical 
GR in the limit when the spacetime curvature becomes small. Quantum evolution in LQC is 
governed by a quantum difference equation, however, for a class of semi-classical states, there 
exists an effective continuum spacetime description. This allows one to obtain an effective 



Hamiltonian constraint for LQC |16Hl8j. The resulting Hamilton's equation lead to the 
modified Einstein's equations which turn out to describe the underlying quantum evolution 
extremely accurately [2|-0, 0]-^ The modified Einstein's equations have been widely used 
in conjunction with the analytical and numerical techniques to capture the details of the 
underlying physics in LQC and to reveal rich phenomenological features fl|. 

In LQC, one starts with a classical phase space in Ashtekar variables: the SU(2) connec- 
tion Al^ and conjugate triads E^, which are then symmetry reduced to c* and pi respectively, 
by incorporating the underlying symmetries of a homogeneous spacetime. The elementary 
variables used for quantization are the holonomies of the connection components and fluxes 
corresponding to the triads.^ The Hamiltonian constraint, the only non-trivial constraint 
left after symmetry reduction, is then expressed in terms of these elementary variables and 
quantized. This procedure leads to two novel features. The first of these arises by express- 
ing field strength of the connection in terms of holonomies over a closed loop shrunk to 
a minimum area on the quantum geometry. This leads to a non-local nature of the field 
strength, which results in a quantum difference equation. The second feature arises due to 
the presence of inverse triad (or inverse volume) operators. Since the eigenvalue spectrum 
of triad operator is discrete and includes zero as an eigenvalue, its inverse it not densely 
defined. Using a classical identity in the phase space, inverse triads are expressed in terms 
of the Poisson brackets between holonomies and the positive powers of triads [20| , and then 
quantized. The resulting eigenvalues of such an inverse triad operator show a significant 
departure from the classical behavior near pi = 0. When the triad component vanish, the 
eigenvalue of such an operator is zero. At larger values compared to the Planck value, it 
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approximates the classical behavior. However, such modifications can only be consistently 
defined when the spatial manifold is compact. In case it is non-compact, the modifications 
to the classical behavior depend on the fiducial volume of the fiducial cell introduced to de- 
fine the symplectic structure. Since this cell is an infra-red regulator introduced to regulate 
infinities occurring due to the underlying non-compactness, physics must be independent 
of it. Indeed, in the limit when the fiducial volume of the fiducial cell becomes infinite for 
the non-compact models, the inverse triad modifications vanish. Thus, for spatial manifolds 
which are non-compact, quantum geometry leads to new physics only via non-local nature of 
the field strength. If the spatial manifold is compact, apart from the non-local field strength, 
inverse triad effects can also be important. 

Let us now discuss some of the main features of the isotropic and Bianchi-I models in 
LQC These are the following: (i) For all of the isotropic and the Bianchi-I models, the 
field strength can be expressed using holonomies, which are almost periodic functions of the 
connection A^, over a closed loop. Quantum geometry fixes the minimum area of the loop 
using the eigenvalues of the area operator in LQG. It results in a universal bound on the 
energy density"^ (Pmax ~ 0.41ppianck) and the expansion scalar in these models, independent 
of the choice of matter (and hence the energy conditions). These bounds have been shown 
to result in a generic resolution of all strong curvature singularities in isotorpic and spatially 
fiat LQC [13] . An analysis of various singularities, including the various exotic ones, strongly 



indicates this result to extendable to spatially curved models j23|. The non-local nature of 



above field strength is also responsible for a universal bound on directional Hubble rates 
and cr^ in the Bianchi-I model, which has been demonstrated to yield resolution of all strong 
singularities for different types of matter [2^; (ii) Inverse triad corrections are potentially 
significant for only those universes which attain a size comparable to Planck value. If bounce 
of the universe occurs when its volume is much bigger than Planck volume, then inverse 
triad corrections play little role. In such cases, they are found to be neither responsible 
for bounds on the energy density and 6, nor do they lead to any significant effects on the 
modified dynamics 0, m. 

This uniformity of results for isotropic models with different spatial curvatures and 
Bianchi-I model is noteworthy. We recall that for the k = 1 model, construction of the 
closed loop is technically challenging where one can not simply extend the strategy used in 
the k = model. ^ The resulting quantization leads to non-trivial terms in the quantum 
Hamiltonian constraint which capture the spatial curvature. In contrast, the quantization 
of A; = — 1 model follows neither the approach in k = model nor k = 1 model to express 
the field strength in terms of holonomies. This is due to a technical difficulty resulting 
from the presence of the off-diagonal terms in the spin connection FJ^. Currently available 
quantizations overcome this problem by considering holonomies of the extrinsic curvature 



Ul, Il2[, and a priori an agreement on the detailed physics at the Planck scale is unexpected. 
Further, in the Bianchi-I model, freezing of the anisotropic degrees of freedom does not 
lead to the quantization of the isotropic flat model. Instead one must integrate out the 
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physical Hilbert space in an exactly soluble model in the spatially flat case 21 1 . 
^ In the spatially flat model, construction of such a loop is straightforward due to the availability of com- 
muting left invariant vector fields on the spatial manifold. In the spatially closed model, left invariant 
vector fields do not commute, and the loop is constructed using both left and right invariant vector fields 
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anisotropic degrees of freedom [13|. Despite various differences in the quantization strategy 
of these spacetimes and the resulting subtleties, it is rather remarkable that loop quantiza- 
tion of isotropic models and Bianchi-I spacetime, reveals the same bounds on the physical 
quantities independent of the energy conditions and are also not affected by inverse triad 
modifications for spacetimes with non-vanishing spatial curvature. 

The goal of this manuscript is to analyze the physics of Bianchi-II and Bianchi-IX space- 
times in LQC, in the above context using the effective Hamiltonian approach. In these 
models, due to the interplay of spatial curvature and anisotropics, it is not possible to 
express field strength in terms of holonomies which are almost periodic functions of the con- 
nection components. A new strategy is needed to loop quantize these spacetimes, which was 



proposed in Ref. [13|. One expresses the field strength in terms of a non-local connection 
defined via holonomies over open segments.^ The underlying quantum geometry does not 
directly constrain the length of such open segments. However, by demanding that the re- 
sulting expression for the field strength operator agrees with the one in the Bianchi-I model, 
one can fix the minimum length of the edge over which the holonomy is computed using the 
minimum quantum of area in LQG. The pertinent question is whether this strategy leads 
to physics at Planck scale which confirms with that established by the quantization of all 
other models in LQC or does it lead to novel surprises? In particular, are energy density, 9 
and cr^ bounded in Bianchi-II and Bianchi-IX models? Are these bounds universal? What 
are the contributions of the modifications originating from the inverse triad operators for 
the Bianchi-IX model? Are these modifications important or are they insignificant, as in 
the case for the isotropic models? 

We will show in this manuscript that the connection operator approach leads to some 
unexpected results. We find that energy density and the expansion and shear scalars are 
bounded in Bianchi-II model, however these bounds are not universal as they depend on the 
energy conditions.^ In the Bianchi-II model, one must assume that the energy density is 
bounded below, else the shear scalar diverges in the Planck regime. In the Bianchi-IX model, 
inverse triad modifications turn out to be critical to obtain these bounds when the approach 
to classical singularity is isotropic, i.e. all triads approach the singularity at the same time. 
For the other types of singularities, we show that energy density and the expansion and 
shear scalars are finite except for two cases. These two cases correspond to the vanishing 
of one or two of the three triads. At this stage, it is neither evident whether the effective 
dynamics can be trusted in such a regime nor if such solutions exist in the Bianchi-IX 
model for some matter content. These results bring out contrasts between the connection 
operator approach [li, 15| and the holonomy approach 0-1, i, 0, H [H, [13 . We find that 



demanding the consistency of connection operator approach with Bianchi-I model, which 
has vanishing spatial curvature, does not guarantee its consistency for the spatially curved 
models. In particular, the isotropic limit of the effective Hamiltonian constraint in the loop 
quantization of Bianchi-IX model does not lead to the effective Hamiltonian constraint of 
the k = 1 model [13], but to that of a different quantization - the inequivalent connection 
operator quantization oi k = 1 model [1^, and that too only in the regime where inverse 
triad modifications can be ignored. Unlike the holonomy based quantization oik = 1 model. 



^ In terms of the symmetry reduced connection, in this approach, the connection operator c is defined as: 

- _ sin(pc) 

' In Bianchi-II model, the energy density has been shown to be bounded for the case of a massless scalar 
field fl^ l . Since the matter content was fixed, the role of energy conditions was not known. 
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where expansion scalar is bounded by a global maxima in the effective spacetime description, 
9 does not have a maximum in the connection operator approach. The same turns out to 
be true for an alternate quantization of A; = — 1 model based on the connection operator 
approach. 

The organization of this paper is as follows. We start with a summary of the effective 
Hamiltonian in LQC for Bianchi-I model and revisit the analysis in Ref. [26[ to show 
that the energy density, and the expansion and the shear scalars are bounded by universal 
values for arbitrary matter content. Here we obtain the correct bound on the shear scalar, 
which was previously estimated incorrectly. In Sec. 3, we repeat the analysis for non- 
compact Bianchi-II model. Here we show that the energy density, the expansion and the 
shear scalar are bounded only if one imposes energy conditions. We show that if one allows 
energy density to be unbounded below, the shear scalar diverges. In Sec. 4, we analyze 
the effective Hamiltonian constraint for the Bianchi-IX model and show that inverse triad 
corrections play an important role for p, 9 and cr^ to be bounded for the isotropic approach 
to singularities. Due to the inverse triad modifications, there is an additional subtlety in 
the behavior of energy density, which is also addressed. Here we also discuss the isotropic 
limit of the effective Hamiltonian in Bianchi-IX model, and find that it does not lead to 



the one for k = 1 loop quantization [lOj- In the absence of inverse triad modifications the 
limit is given by an alternative connection operator based quantization, which has expansion 
scalar unbounded. Whereas in the presence of inverse triad modifications, the isotropic limit 
agrees with none of the available isotropic quantizations. We summarize the results with a 
discussion in Sec. 5. In Appendix A, we discuss the behavior of expansion scalar in = 1 
model in the holonomy based and the connection operator quantizations in LQC. A similar 
discussion for k = —1 model is provided in Appendix B. These show that for isotropic 
k = ±1 models, the expansion scalar is unbounded for the effective Hamiltonian constraint 
corresponding to the connection operator approach. 



II. BIANCHI-I MODEL 

Bianchi-I model is one of the simplest examples of an anisotropic spacetime. It has van- 
ishing spatial curvature and in the isotropic limit it yields the k = Robertson- Walker (RW) 
metric. The homogeneous Bianchi-I anisotropic spacetime can be described by a manifold 
S X M where E is the 3-spatial hypersurface characterized by a set of three commuting Killing 
vectors. As earlier works in LQC, we will consider the topology of S as M^. The spacetime 
metric of Bianchi-I spacetime is given by, 

ds"^ = -N'^dr'^ + a\dx^ + aldy'^ + aldz^ (2.1) 

where ai,a2 and are the directional scale factors. These can be used to define a mean scale 
factor a := (010203)^/'^. Since the spatial manifold is non-compact, in order to define the 
symplectic structure, one needs to introduce a fiducial cell V. This cell's edges can be chosen 
to lie along the integral curves of fiducial triads and have coordinate length li,l2,h- The 
cell V has fiducial volume Vo = hhh with respect to the fiducial metric qab compatible with 
fiducial co-triads w*. 

In LQG the canonical variables are the Ashtekar connection and the triad Ef which, 
which due to the symmetry of Bianchi-I spacetime, can be expressed in terms of components 
c* and Pi as: 

Ai = d{h)-'ul and = PikV-^^q (2.2) 
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The connection and triad components, c* and pi, satisfy the following Poisson bracket rela- 
tion, 

{c\p,} = 87cG-f6] (2.3) 

where 7 ~ 0.2375 is Barbero-Immirzi parameter. The triads pi are related to the directional 
scale factors as 

Pi = 61^2^3 |«2a3|, P2 = £2^1 ^3 |ai osl, Pa = £3^2^i|aia2| (2.4) 

where Si — ±1 depending on the orientation of the triads. Without any loss of generality, 
we will choose the orientation to be positive in the following analysis. 

Let us first consider the Hamiltonian constraint in the classical theory. For lapse N = V, 
the classical Hamiltonian constraint in terms of Ashtekar variables can be written as. 



"^ci = o \, o (C1P1C2P2 + cychc terms) + TimattV , (2.5) 



where "Hmatt is the matter Hamiltonian and V denotes the physical volume of the cell V: 
V = v'PiP2P3- Dynamical equations can be obtained using Hamilton's equations: 

Pi = {p^,n,l} = (2.6) 

and 

Ci = {ci,n,i} = 87rG7^ (2.7) 

OPi 

where the 'dot' represents the derivative with respect to the proper time t. Using the first 
equation one can compute the directional Hubble rates, Hj — ai/ai, such as 

2 \P2 P3 Pi 

and show that the connection component q are related to Hi in the classical theory as 
Q = jliHiai. These directional Hubble rates also define the expansion scalar 9 and the 
shear scalar for the comoving observers: 

^= 1^ = (H^ + H^ + Hs) (2.8) 

and 

a' = J2 (H^ -^f {{Hi - H^f + {H2 - H,y + {Hs - H,f) . (2.9) 

i=l 

In an isotropic spacetime Hubble rates in all directions are equal and the shear scalar van- 
ishes. 

Physical solutions of the Einstein's equations satisfy the constraint l-Ld ~ 0, which results 
in the following equation for the classical GR: 

HrH2 + H2H^ + = SttG p (2.10) 

where the energy density p is computed as p = Hmaxt/y ■ Using the equations for i^j's, it is 
possible to obtain the equation for the mean Hubble rate H — a/a, which turns out to be of 
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a similar form as the Friedmann equation in the flat isotropic model albeit with a presence 
of a terms proportional to the anisotropic shear: 



"'—"^^ (2.11) 

where T? := |(j^a^. In the classical theory, turns out to be a constant of motion. Note 
that due to the presence of positive deflnite shear scalar, the energy density can take negative 
values without Hubble rate becoming imaginary. This is in contrast to the isotropic limit of 
the above equation, corresponding to the /c = FRW model, where negative energy densities 
are not allowed when Hubble rate is real. Using above dynamical equations, one flnds that 
at vanishing scale factors, p, 9 and diverge. These lead to the divergence in curvature 
invariants and the break down of geodesic evolution at the singularities.^ 

Let us now discuss the dynamics in LQC based on the effective Hamiltonian approach 
based on the embedding method [igI- 18|.^ This approach is based on geometrical formulation 



of quantum mechanics [29}], where the space of quantum states is regarded as an infinite 
dimensional quantum phase space, Fg, equipped with a symplectic structure provided by 
the Hermitian inner product on the Hilbert space. Using a judicious choice of states, one then 
finds an embedding of the finite dimensional classical phase space F. If this embedding is 
approximately preserved under the fiow generated by the quantum Hamiltonian vector field, 
then up to the order of approximation the embedding is considered faithful, and one obtains 
effective equations which incorporate quantum corrections. In LQC, this procedure has been 



successfully carried out for a variety of models, including a dust dominated universe [16 



isotropic model with a massless scalar field [17| and isotropic model with arbitrary matter 



[18[. The resulting effective Hamiltonian leads to modified Einstein's equations which have 
been extensively compared with the true quantum evolution using numerical simulations 



with states which are semi-classical at late times [3H6|, [lOl, lUl • These analyses show that 
the modified Einstein's equations obtained from the effective Hamiltonian in LQC are in 
excellent agreement with the quantum evolution for all the models considered so far. Thus, 
suggesting that the effective Hamiltonian approach in LQC may have validity for a large 
class of spacetimes. 

The effective Hamiltonian constraint for Bianchi-I spacetime in LQC for lapse = is 



given as [13|, [261, |30 



1 /sin (uiCi) sin (U2C2) \ ^, , „x 

^cfi = ^^PiP2 + cyclic terms + V,^^,,V , (2.12) 



To give an example, we recall that the square of the Weyl curvature C'a/3^1, can be expressed in terms of its 
electric {Eap) and magnetic {Hap) parts relative to a unit time-like vector field (u): Eap = Cc^ypsu^u^ , 
Hcfi = *Ca^t0su<u\ as Ca^^^C"''^'^ = ^{EapE'^^ ~ H^f^H'^P) (see for eg. 0). For the Bianchi-I model, 
Elap — and Eap — + fcQ/s — (c^c-y^ — ^a^Sap) when anisotropic stress vanishes. Thus, a 

divergence in expansion and shear scalars, leada^to a divergence in the square of the Weyl curvature. 
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Another approach is the truncation method [28[, also based on geometric formulation of quantum me- 
chanics. In contrast to the embedding approach which can be regarded as analogous to the variational 
methods, the truncation approach is on the lines of order by order perturbation theory. Though the 
method is more systematic, one needs to exercise a lot of care in dealing with truncation errors. Un- 
like the embedding approach, the modified equations from this method have not been widely tested for 
compatibility with the underlying quantum evolution. 
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where /i, are given by 



/ii = A./-^; /i2 = A./-^;and jj,^ = X (2.13) 
V P2P3 V P^P3 V 

and A^ = Ay/STi'jlpi. Note that "Hmatt here corresponds to the matter Hamihonian as obtained 
from the Fock quantization. The relationship of fli with the triads is a consequence of the 
way area of the loops along which holonomies of connection are computed relates with the 
eigenvalue of the area operator in LQG.^° The minimum allowed area of such a loop is 
labeled by A^. It is to be noted that above functional dependence of /i, on the triads is 
unique, in the sense that any other choice leads to resulting physics being affected by the 
rescaling of the lengths of the edges of the fiducial cell and also by change in its shape 



31[. Further, the only contribution from quantum geometry in the effective Hamiltonian 
constraint results from the modifications originating from the non-local field strength. The 
inverse volume modifications are absent in the above constraint, since we are considering a 
non-compact spatial manifold. 

Using Hamilton's equations, we can compute the time variation of triads. In terms of 
proper time, we get 

Pi - - - 

Pi = — (sin (/i2C2) + sin (fisCs)) cos (fiiCi) (2.14) 
7A 

and similarly for p2 and ps, and the connection components. From these, we can obtain the 
equations for the directional Hubble rates, such as 

Hi = (sin {fiici - 112C2) + sin (/i2C2 + /iaCs) + sin {fiiCi - flsCs)) (2.15) 
27A 

Using these equations, one finds that that unlike the classical theory, q 7^ 'jljHjai. Further, 
it is straightforward to show, that is not a constant of motion in LQC [261, |30 . 



In contrast to the classical theory, the directional Hubble rates are bounded in LQC, 
with a maximum value given by -ffj^max = 3/27A. The resulting expression for the expansion 
scalar yields 

= (sin (/iici + /i2C2) + sin (/i2C2 + fisCs) + sin {fiiCi + fisCs)) (2.16) 
27A 



which has the following maxima 



27A /pl 



The boundedness of the Hubble rates points to the non-singular bounces in LQC. These 
bounces, unlike the isotropic case, do not occur at a fixed value of energy density. Due 
to the presence of anisotropies, bounce occurs at different values of p and cr^, which are 
determined by the initial conditions (see for example, Ref. jlHl)- In order to find the 



It is important to note that at the current stage of research, the imposition of the minimum area of the 
loop as it appears in LQG is an external input in LQC. Since a derivation of LQC from LQC is yet to 
be performed, LQC is a quantization of cosmological models based on LQG, rather than the cosmological 
sector of LQG. 
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maximum values, we first note that the vanishing of the Hamiltonian constraint, HcS ~ 0, 
gives 

P = (sin (/iiCi) sin (/iaCs) + cyclic terms) , (2.18) 

which implies that p has a universal maxima, independent of any energy conditions, as: 

Pn.ax = ^^^J^ ^ 0.41PP1 . (2.19) 

Thus, the maxima of energy density in Bianchi-1 model turns out to be the same as in the 
isotropic models in LQC Note that for an arbitrary choice of initial conditions in Bianchi-I 
spacetime, is never achieved at a bounce of the directional scale factor. This bound is only 
saturated when anisotropies vanish. 

Similarly, using eq. fl2.15p (and equations for H2 and if 3) in fl2.9p . it is straightforward to 
obtain the expression for shear scalar, which as the one for 9 and p turns out to be composed 
of bounded functions: 

or = „ (cos (/i2C2)(sin (pici) + sin (psCs)) - cos (/iiCi)(sin {P2C2) + sin (/^aCs)))^ 



' 372A2 

+ cyclic terms 



(2.20) 



(where the subscript / is used to differentiate this expression from shear scalars in subsequent 
models considered here). The shear scalar has a global maxima 



_ 10.125 ^ 11.57 



at piCi = 7r/6,/i2C2 = 7r/2 and p^c^ = 57r/6. Interestingly, the maxima of shear scalar is 
reached when energy density is itself close to the Planckian value. At above values of PiCi, 
the energy density turns out to be p ~ 0.4167pcrit- In the case of Bianchi-1 vacuum spacetime 
(p = 0), depending on the initial conditions, the shear scalar can attain one of the two (local) 
maximas: cr'^r^^ = 2/7^A^ at piCi = tt/2,P2C2 = vr and P3C3 = 0, or cr^g^ = 6.3/37^A^ at 
PiCi = —0.339837, P2C2 = 7r/2 and P3C3 = 5n/6. 

The existence of upper bounds on energy density, expansion scalar and shear scalar 
for arbitrary matter indicate a resolution of various singularities in the loop quantization of 
Bianchi-1 spacetime. This is a direct feature of the underlying quantum geometry captured 
via A^. In the limit, this parameter goes to zero, one recovers the classical divergence of p, 9 
and (T^. this feature is missing in the classical theory. The upper bounds of these quantities 
are completely generic, obeyed by all types of matter. As we will discuss in the next sections, 
this generality is absent in Bianchi-11 and Bianchi-lX models. 



Strictly speaking, a stronger statement can be made for the isotropic models in LQC. This is due to 
the reason that in the isotropic case, availability of an exactly soluble model [2l|, allows one to prove 
existence of bound on energy density for a dense set of states in the physical Hilbert space. Where as 
here, bounds are derived assuming an effective Hamiltonian which a priori assumes coherent states. Thus, 
it is possible that above bounds, and similarly those derived later for Bianchi-II and Bianchi-IX models, 
are not necessarily strict bounds for all the states in the physical Hilbert space. 
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III. BIANCHI-II MODEL 



Bianchi-II spacetimes are more general than the Bianchi-I spacetimes in the sense that 
they have a non-vanishing intrinsic curvature. Unhke the Bianchi-I model, only two of the 
three Killing fields commute with each other on the spatial manifold. Further, the spacetime 
lacks an isotropic limit. In the following we consider the case of a non-compact Bianchi-II 
spacetime. Since the procedure of symmetry reduction, introduction of a cell V and details of 
the symplectic structure follow closely as in the Bianchi-I model, we will skip the discussion 
of this part. For a detailed discussion on these for Bianchi-II model, we refer the reader to 
Refs. Q. 

In terms of Ashtekar variables, the classical Hamiltonian constraint for lapse N = V, can 
be written as 



^cl 



1 



[P1P2C1C2 + cyclic terms] 



1 



ap2P3Ci - (1 + 7 ) I ) 



(3.1) 

where a is related to the structure constants of the Lie algebra corresponding to the Killing 
fields. Using Hamilton's equations, it is straightforward to obtain the classical dynamical 
equations for the time variation of Pi and Cj. As an example: 



dpi 
dr 



1 
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(Pl(c2P2 + C3P3) + aP2P3) 



(3.2) 



However, unlike the Bianchi-I model, it is not possible to use these equations to obtain an 
analog of the generalized Friedmann equation (12. lip for the mean Hubble rate. A detailed 
analysis of these dynamical equations reveal a singularity when the triads vanish. 

We now consider the effective Hamiltonian for the Bianchi-II model [l3]. For the same 
choice of the lapse function as above, it is given as 



H 



eff 



P^P'^Pj^ [sin (^^ci) sin {H2C2) + cychc terms] 



a{P2P3f''^ . I- ^ 

-sm (/iici, 



^2(1+72 



A 



P\ 



+ nm.ttV (3.3) 



which is a sum of the effective Hamiltonian in Bianchi-I model and the terms originating 
from the presence of the spatial curvature. Imposing the the constraint Hcs ~ we obtain 



1 



P 



[sin (/iiCi) sin (/i2C2) + cyclic terms] 



where as in the Bianchi-I model, p = 'Hma.tt/V and x is defined as 



X . , . (1 + 72)^2 
- sm (/iici) ^ 



X = a 



P2P3 
Pi 



(3.4) 
(3.5) 



12 The Killing fields satisfy [^i.^j] = C^j^k, with the only non-zero structure constant being C23 
This defines a as a := {l2l3/h)ct, where k refer to the edge lengths of the fiducial cell V. 



a. 
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The right hand side of the eq. (13.40 has a global maxima at x 
This results in a maximum allowed value of the energy density as [13 



and sin(/ijCj 



1. 



3 + (l + 7^) 



OMp 



PI • 



(3.6) 



The above upper bound on p is the global maxima. It is straightforward to see that eq. 
(13. 4p can give rise to large negative energy densities depending on the values of triads. 
This behavior is depicted in Fig. 1. A possible way to obtain a lower bound on the energy 
density is by imposing energy conditions on the matter content. As an example, if we assume 
weak energy condition (WEC), then p > 0. This procedure is similar to imposing energy 
conditions in GR to eliminate solutions with negative energy density. However, unlike GR, 
where even after imposing energy conditions, energy density is not bounded above, in LQC 
imposition of an energy condition on p does not affect the global maxima of p which as we 
have shown above arises because of the underlying quantum geometric effects. 




FIG. 1: The plot (in Planck units) shows the variation of energy density as a function of x in the 
loop quantization of non-compact Bianchi-II model when sin(/ijCj) = 1 is substituted in eq. (j3.4p . 



Using (13. 3p . one can derive the modified dynamical equations using Hamilton's equations. 
For the triad components we obtain 

dpi 1 f jy^ \ 

— = - — (sin (/isca) + sin [p^c-i]) + apaPs cos {piCi) 

dr 7 V^i / 

which yields the classical equation (13. 2 p in the limit — )■ 0. In terms of the proper time, 
this equation can be written as 



Pi \ _ _ _ 

— = — (sin {P2C2) + sin (/isCs) + Ax) cos {piCi) . 

Pi iX 



(3.7) 



As noted below eq. (|2.1ip . in the anisotropic spacetimes, negative values of energy density do not necessarily 
cause problems with the reality of the expansion scalar. 
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Equations for the time variation of other triads can be derived in a similar way, and they 
turn out to be: 



P2 1 

— = — (sin (/2iCi) + sin (/23C3)) cos (/iaCs) (3.8) 
Ps 1 

— = — (sin (/iici) + sin {^202)) cos {ftsCs) ■ (3.9) 
Ps 7^ 

Using these equation we can obtain the directional Hubble rates, Hi, and the expansion 
scalar 6 using eq.f l2.8p which becomes 

I _ _ 

9 = —— (sin(/iici + /X2C2) + sin(/i2C2 + /xsCs) + sin(/X3C3 + /xiCi) + Ax cos (/xiCi)) . (3.10) 

Unlike the expansion scalar in the Bianchi-I model, 6 in Bianchi-II model is generically 
unbounded because of the divergence in a; as pi — )■ or {p2,P3 00). However, if one 
imposes energy conditions demanding that the energy density be bounded below, then x can 
not grow beyond a maximum value, and 6 turns out to be bounded. For matter satisfying 
WEC, the maximum allowed value of 6 is given by, 

,3.n) 

27A Ipi 

This upper bound can be obtained by finding the maxima of eq.f l3.10p numerically by op- 
timizing various variables and occurs at fliCi = 0.642, p,2C2 = 0.982, /isCa = 0.982 and 
X = 1.717. 

Finally, the expression for the shear scalar in Bianchi-II model, Uj^j can be obtained using 
(12. 9p . which after a straightforward calculation yields. 



^2 ^2,1 

CTu = O-j + 



2\^x^ cos^(/iiCi) + 2Aa; cos (/iiCi)(2(sin (/X2C2) + sin (/xsCa)) cos (fiiCi) 

(3.12) 



II - -I 

-(sin (/iici) + sin (fxsCs)) cos (/22C2) - (sin (1x202) + sin (fxiCi)) cos (/iaCs)) 



where ai is given by eq.f l2.20p . This term is bounded, as shown in Sec. II. However, the term 
with square parenthesis is an increasing quadratic function in x. Thus, cx^j does not have 
an upper bound. The shear scalar for non-compact Bianchi-II model is devoid of a generic 
maxima in contrast to the Bianchi-I spacetime. In order to obtain a bound on (jjj, one needs 
to assume a lower bound for the energy density, as for the expansion scalar. Imposing WEC, 
the maximum value of cxjj turns out to be 



2 57.58 65.82 

^11 max ^ ^ ^ ■ (3.13) 



Note that the maximum of the shear scalar, as that of the energy density and the expansion 
scalar in the non-compact Bianchi-II model turns out to be different from the Bianchi-I 
model. 



12 



It is to be emphasized that the bound on the shear scalar is sensitive to the energy con- 
ditions, i.e. if one imposes a different energy condition, the allowed range for the parameter 
X would be different and so would be the upper bound on the shear scalar. WEC is followed 
by almost all types of the matter and it provides an adequate estimation of the bounds on 
shear and energy density. Thus, WEC has two roles: first to give a lower bound on p and 
the second to provide the corresponding upper bounds on 9 and (Jjj. Further, the upper 
bounds on 9 and cr^ are saturated in the dynamical evolution, and the bounce of directional 
scale factors occur before these values are reached. As in the case of Bianchi-I model, this 
is tied to the interplay of the Ricci and the Weyl parts of the spacetime curvature in the 
Bianchi-II model. 



IV. BIANCHI-IX MODEL 

The spatial manifold for the Bianchi-IX model has a compact topology whose radius 
with respect to the fiducial metric is chosen as = 2. The volume of fiducial cell is given 
by Vo —: il — 27r^a^ = 167r^. In contrast to the Bianchi-I and Bianchi-II spacetimes, none 
of the Killing vectors commute with each other in the Bianchi-IX model. In this spacetime, 
the interplay of intrinsic curvature and anisotropics is also much richer in comparison to the 
Bianchi-II model. Classical dynamics can exhibit a Mixmaster behavior as singularities are 
approached. In the isotropic limit, one recovers the classical dynamics of the k — 1 FRW 
model. 

Let us first consider the classical Hamiltonian constraint for the Bianchi-IX model. With 
lapse N = V, in terms of the connections and triads, it is given by 



f)i 



2£o 

P1P2C1C2 + P2P3C2C3 + P3P1C3C1 + -Y (P1P2C3 + P2P3C1 + P3P1C2) 



P3PI 
P2 




which using Hamilton's equation for motion, leads to 



dpi 
dr 



Pi 

7 



P2C2 + P3C3 + io 



P2P3 
Pi 



(4.2) 



and similarly for p2 and ps, and the connection components. In this case, one can use these 
equations to derive a generalized Priedmann equation, as in the Bianchi-I model, and it 
turns out to be: 



= p+ -a^ 

3 6 12pip2P3 



2{pI+pI+pI)- 



P1P2 

P3 



P2P3 
Pi 



P3P1 
P2 



(4.3) 



The Hubble rate diverges as the singularities are approached in the Bianchi-IX spacetime. 
Further, due to the presence of intrinsic curvature, singularities can occur both in past and 
the future evolution, as in the k = 1 universe. As shown in Appendix A, in the isotropic 
limit, one recovers the dynamical equations for the k — 1 model. 
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As in the case of the loop quantization of the Bianchi-II model, it is not possible to express 
the field strength operator in terms of holonomies along closed loops. To overcome this dif- 
ficulty, a quantization has been proposed, along the lines of Bianchi-II model, following the 
connection operator approach [l5|. In the following we analyze the physics resulting from 
the effective Hamiltonian constraint of this quantization. Since the underlying manifold 
is spatially compact, the resulting effective Hamiltonian constraint contains modifications 
originating from both the non-local nature of the field strength operator and the eigenvalues 
of the inverse triads. Let us first analyze some features of the resulting physics for the effec- 
tive Hamiltonian constraint of Bianchi-IX model, if inverse triad modifications are ignored 
(as in the analysis of Ref. [lH). We denote this effective Hamiltonian constraint by Hcs in 
order to distinguish it from the Hamiltonian constraint in eq. (14. lip where the inverse triad 
corrections are included. For the lapse N = V, the effective Hamiltonian constraint is given 
by Q: 



^eflf = o"^!^^^^ o [sin (yWiCi) sin {^1202) + cychc terms] - 



327rG72 



3/2 



2{p\ + p\+ p\) — { ^ ^^^^ ^ + cyclic terms 



■ sin (/iaCs) + cychc terms 



+ V-r..ttV . (4.4) 



The vanishing of this constraint leads to the following expression of energy density 
1 



P 



o ^ TT? N'^ {PiCij sin {H2C2) + cyclic terms] + 



^ sin (;Z3C3) + cyclic terms 



327rG72 



,Pi_ 

P2P2. 



P2Pz \ 
Pi ) 



cyclic terms 



(4.5) 



Utilizing the boundedness properties of the trigonometric functions, this expression implies 
that 



+ 



87rG72A2 87rG'72A 
where xi = 



[Xi +X2 + a^3] + ^°o„^„2 [2(XiX2 + XaXg + X3X1) - (x^ + + xf)] 



3271^72 



(4.6) 



P2P3 



X2 



PlP3 



and X3 



_ I P1P2 

Pi ' ' \ P2 ' """" V 'W' 

The first term in eq. (14. 6 p is the maximum energy density obtained in Bianchi-I spacetime. 
It is the behavior of the second and the third term, which determines the boundedness of 
energy density from the above effective Hamiltonian constraint. In order that the energy 
density to have a maxima, p{xi,X2,X3) should have a viable simultaneous solution to the 
following system of equations: 



dp{Xi,X2,Xs) 



0; 



(9p(Xi,X2,X3) 



^ dp{Xi,X2,X3) ^ Q 



dxi 8x2 9x3 

Solving these equations, one obtains the following condition: 

6 



^oA(l + 7' 



+ Xi + X2 + Xs = 



(4.7) 



(4.8) 
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Thus a physical solution is allowed only when at least one of Xi is negative. However, 
by definition, all of Xi are positive}'^ The energy density resulting from Hes does not 
have a maxima, unlike the Bianchi-I and Bianchi-II models in the absence of inverse triad 
modifications. On analyzing the modified dynamical equations resulting from "Heff, we find 
that the expansion and the shear scalars are also unbounded. 

However, as discussed above the effective Hamiltonian constraint "Hefr is incomplete as 
it lacks the contribution from inverse triad corrections in the Bianchi-IX model which is 
spatially compact. To conclude whether the energy density in Bianchi-IX model has an 
upper bound, it is necessary to include these modifications. To introduce these modifications 
in the corresponding effective Hamiltonian constraint in LQC for such terms, we consider 

eigenvalues f{pi) of the operator Pi~^^'^ and substitute them in place of inverse triad terms. 

As an example, the eigenvalues for the inverse triad operator pi~^/^ turn out to be [JJ 



Pi \Pl,P2,P3) = fiPl)\pi,P2,P3) 

where \pi,P2,P3) denote eigenstates of the volume operator, and v 

2 



(4.9) 

^ T-y/PiPm and 



fiPi 



47r7A/pi 



{P2P3) 



1/2 



f + 1 



V - 1 



(4.10) 



Similarly, one can derive expressions for other inverse triad operators. With these modifica- 
tions, the effective Hamiltonian constraint becomes: 



n 



cff 



P1P2P3 

eio 

1(1+ 7') 



sin (/iiCi) sin (/i2C2) + cychc terms 



(PlP2)^^V(P3) sin (/igCs) + {P2P3T'^ fi.Pl) sin (/iiCi) + {p3PlT^^fiP2) sin (/i2C2 



\3/2 



3/2 . 



327rG72 



2{pI+pI+pI) - {pm?f{p3r - {P2P3?f{PlY - {P3Plff{P2) 



(4.11) 



where the inverse triad modifications also contribute to "Hmatt if one considers matter Hamil- 
tonian containing inverse powers of the scale factor. 

We now analyze the behavior of energy density. Due to the inverse triad modifications, an 
ambiguity in its definition arises. Using the effective Hamiltonian constraint, one can define 
the energy density p as p = Timatt/V, where the "Hmatt includes modifications due to inverse 
triad operators. However, one can also define energy density such that it agrees with the 



The unboundedness argument given here is vahd irrespective of the choice of orientation of triads, which 
has been here fixed to be positive. 

The strategy is same as in the isotropic models and is based on the identities first proved by Thiemann in 
LQG [2^ . For an up to date discussion of these modifications in isotropic models and related subtleties 
for non-compact models, see Ref. (if. 



15 



eigenvalues of V'^T-Lmnu (suitably symmetrized). We label this energy density as Pg. In the 
absence of inverse volume modifications, p and Pg are equal to each other. However, when 
these are present, Pq and p can behave in a qualitatively different way, at small volumes. 
As it turns out, p and Pq, indeed have qualitative differences in the loop quantization of 
Bianchi-IX model. 

Using the effective Hamiltonian, a division by volume, leads to the following inequality 
for the energy density p for the physical solutions: 



P 



< 



+ 



87rG72A2 87rG72A 



327rG72 



Pi 



y/PlP2 



+ 



P3 
P2 



Pi 



+ 



P3 



P2P3 P1P3 P2P1 



P1P2 r( X4 

p-i 



P2 

_ P2P3 
Pi 



(4.12) 



fiPl 



PsPi 

P2 



f{p2r 



p 

50- 



FIG. 2: Behavior of energy density p as a function of v is shown (in Planck units). 



The behavior of energy density with respect to v is shown in Fig. 2. We find that the 
energy density p does not have a global maxima if the effective Hamiltonian description 
is assumed to be valid for the entire range of v. However, in the case when all the triads 
approach the classical singularity at the same time i.e. the singularity is isotropic, then for 
the range f > 1 there exists a local maxima at pi ^ 3.634/pj , given by 



Pn 



6.34ppi . 



(4.13) 



Thus, for the energy density defined as the ratio of the matter Hamiltonian to the physical 
volume, inverse triad modifications do not suffice to control the divergence as f — )■ 0. 

Let us now analyze the behavior of pq. In this case, the effective Hamiltonian constraint 
yields 



■^^ In LQC, physics has been analyzed using both of these definitions. For a comparison of some of the 
features, see Ref. 32 1. 
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87rG72A2 
■^(1 + 7^) 



^7rG72A 



P1P2 r, X , P2P3 rt X , Pm ft 



327rG72 
{j)2Pzf''^ 



3/2 
Pi 



+ 



P3 

3/2 
P2 



+ 



3/2 
P3 



Pi 

{PlP2f^^ 



(4.14) 



fipsY 



fipiT 



where g{v) denotes the eigenvalue of the volume operator (l/V^) 

1 



\/v + 1 — y/v — 1 



(4.15) 



We find that there exists a global maxima at pi ^ 2.109/pi, when the approach to classical 
singularity is isotropic, given by 



11.74ppi . 



(4.16) 



The variation of pq with respect to v is shown in Fig. 3. Note that this maximum value is 
higher than the one in Bianchi-I or Bianchi-II models. When the approach to the classical 
singularities is not isotropic, such as for pancake or cigar singularities, above bounds can be 
violated. This is evident from Fig. 4, where we see that if in the effective dynamics two of 
the triads vanish and the third diverges, then the energy density Pg diverges. In all other 
cases, it is finite. 



4 6 
V 



FIG. 3: Variation of energy density pq is shown versus v in Planck units when the approach to 
classical singularities is isotropic. The minimum occurs at — 753.09ppi. 

Let us now obtain the dynamical equations for the triads and the expression for the expan- 
sion scalar. Using the effective Hamiltonian constraint (14. lip , one obtains using Hamilton's 
equations: 



Pi 
Pi 



1 

7A 



sin {P2C2) + sin (psCs) + e Xio^^^^^fiPi] 

Pi 



cos {piCi] 



(4.17) 
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and similarly for p2 and J93. Using these, the expansion scalar can be computed as 



27A 



sin (/i2C2 + /^scs) +eUoCos {liiCi)^^^^^ f [pi] 

Pi 



cyclic terms 



(4.18) 



The expansion scalar turns out to have a global maxima at Pi 
approach to classical singularity, given by 



2.258/|j, for the isotropic 



6 ^ 



47.72 44.18 



27A 



PI 



(4.19) 



whereas for more general approaches to singularity such as cigar or pancake like ones, there 
exists no such global maxima for 9. 

Note that unlike the Bianchi-II model, this is a global bound which does not depend 
on the imposition of any energy condition. Nevertheless, if one imposes WEC as in the 
Bianchi-II case, one get a lower value. The maximum allowed value of 9 for matter satisfying 
WEC (in terms of Pq) is ^^ax = 30.84/27A. 




FIG. 4: Variation of energy density pq, in Planck units, is shown versus p2 and (with pi = P2). 
The energy density diverges as two of the triads tend to zero while the other approaches infinity. 
For other trajectories pq is finite. 

We now consider the behavior of shear scalar in the Bianchi-IX model. Using Hamilton's 
equations, we obtain 



(T- 



IX 



372A2 



sin {H2C2) + sin (/igCg) + e Xio ^^^^^ fiPi) ] cos {piCi] 

Pi 



VPiP3 



sin (/iiCi) + sin (/^sCs) + A^q-^^ /(P2) ) cos {^202) ) + cyclic terms 

P2 



(4.20) 



For the isotropic approach to classical singularity, the shear scalar has a global maxima at 
Pi ^ 2.258/pi, given by 

2 2165.91 2476.04 

- ~ . (4.21) 



"^IXmax Q_2 \2 



372A2 
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For matter satisfying the WEC (for p,), the maximum allowed value is given by 0"^^^ ~ 
690.98/(37^A^). If the approach to the classical singularity is not isotropic, then the shear 
scalar can take values greater than the above value. This behavior can be seen from Fig. 
5. There exist two cases where shear scalar can diverge in this model, //effective dynamics 
allows solutions where one triad vanishes and two diverge or where two of the triads vanish 
and the third one diverges, then the shear scalar can diverge. It is important to note that the 
existence of above divergent cases is based on the assumption that the effective Hamiltonian 
(14. lip remains valid as above anisotropic singularities are approached. Since the validity of 
(14. lip has not been tested with the underlying quantum evolution, it is possible that these 
cases originate in the regime where the effective dynamics resulting from (14. lip may break 
down. 



10 P3 




FIG. 5: Variation of shear scalar cj^ for Bianchi-IX is shown versus p2 and ps (with pi = P2)- For 
isotropic approach to singularity there exists a local maxima. 



Thus, we find that in the Bianchi-IX model, energy density (pg), expansion scalar and 
the shear scalar are all bounded for isotropic approach to singularity, if we include inverse 
triad modifications to the effective Hamiltonian constraint. However, energy density (p) 
defined by taking a ratio of the matter Hamiltonian to the physical volume does not have 
a global maxima. It is to be emphasized that in the dynamical evolution these bounds are 
not saturated and in certain situations, bounce of scale factors can occur close to the values 
in Bianchi-I and Bianchi-II models. 



We conclude this section, with a discussion of the isotropic limit of the Bianchi-IX model. 
As mentioned earlier, in the classical theory, the limit is given by the classical Hamiltonian 
constraint corresponding to the k = 1 model. Considering the isotropic limit of classical 
Hamiltonian constraint (14.11) . one obtains 



(iso) 



'^1 



(4.22) 



This corresponds to the classical Hamiltonian constraint of = 1 model (eq. (lAip ) with 
e = —1. The pertinent question is whether the effective Hamiltonian constraint in Bianchi- 
IX model leads to the effective Hamiltonian constraint of the k = 1 model in LQC The 
quantization oik = 1 model, based on expressing field strength in terms of holonomies over 
closed loops was performed in Ref. 10|. A detailed analysis of the physics of this model 



reveals features similar to the /c = model in LQC, with a bounce occurring at p = pn 
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0.41ppianck, and resolution of various strong curvature singularities [23|. In Appendix A, we 
discuss the effective Hamiltonian constraint corresponding to this quantization (eq. (1A6P ). 
and show that the expansion scalar in this quantization is bounded above by a universal 
value. 

Let us first consider the effective Hamiltonian constraint without the inverse triad mod- 
ifications. The isotropic limit of "Hefi (eq.f l4.4p ). obtained by imposing pi = p and Cj = c, 
leads to 



(iso) 
eflf 



3p2 



This equation can be written as 

3p' 



sin2(/ic) ^ ^ sin(/xc) ^ ^g(l + 7^ 



ft' 



(iso) 
eff 



877^72 



sin (/ic) 



/i 



+ 7 



+ 



(4.23) 



(4.24) 



which does not agree with (IA6p .^^ Hence, the effective Hamiltonian constraint "Hes in the 
isotropic limit does not yield the effective Hamiltonian constraint of the k = 1 model in the 
quantization based on holonomies over closed loops [lo']. 

It turns out that the isotropic limit of fl4.4p corresponds to the effective Hamiltonian 
constraint of an alternate quantization of A; = 1 model based on the connection operator 
approach (eq. (IA10|) ) discussed in Appendix A. This is straightforward to see using f l4.22p 
which under: c — )■ sin(/2c)//i, yields (14. 24 p . However, as discussed in Appendix A, this 
quantization has the following drawback in comparison to the one in Ref. |lo|. If one 
assumes, the validity of the effective Hamiltonian constraint for all f , the expansion scalar 
turns out to be unbounded even after the inclusion of inverse triad modifications. 



Finally, we consider the isotropic limit of the constraint (14. lip . Imposing Pi = p and 
Cj = c, where p and c refer to the isotropic triad and connection variables, we obtain 

3p3 



n 



(iso) 
eff 



87rG72A2 

3277^72 



[sin^ (/ic)] 



3^0 



87rG72A 

[2p2_p2j(p)4] ^-H^^^^V 



[p^fip) sin(/ic)] 



(4.25) 



nor the alter- 



which agrees neither with the effective Hamiltonian constraint of Ref. 
nate quantization in the presence of inverse triad modifications. In summary, the effective 
Hamiltonian constraint of Bianchi-IX model in the isotropic limit does not yield the effective 
Hamiltonian constraint of the loop quantization of = 1 model performed using holonomies 
over closed loops [lo|. Its isotropic limi t, ig noring inverse triad corrections, agrees with the 
alternate quantization of /c = 1 model 25|, in which the expansion scalar does not have a 



maxima. 



V. DISCUSSION 

In this work, we investigate the behavior of energy density, the expansion and the shear 
scalars in Bianchi-II and Bianchi-IX models using the effective Hamiltonian approach in 



It is straightforward to check that this is true irrespective of the choice of orientation of the triads. 
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LQC These spacetimes have a non-vanishing spatial curvature which leads to technical 
difficulties to carry out the loop quantization as in the Bianchi-I model. In particular, it 
is not possible to express field strength in terms of holonomies around closed loops such 
that the holonomies are almost periodic functions of connection components. The resulting 
holonomy operators are not well defined on the kinematical Hilbert space. To overcome 



this problem, a strategy was proposed in Ref. which expressed field strength directly 
in terms of a non-local connection operator, defined via holonomies computed along open 
segments. In the connection operator approach, coefficients of the expression of field strength 



are fixed by demanding consistency with the Bianchi-I model [IJ]. The resulting quantum 
Hamiltonian constraint turns out to be non-singular as in the Bianchi-I model. However, 
physics of these spacetimes was largely unexplored. First, it was not known whether the 
energy density and the expansion and shear scalars are bounded in these models. It was 
also not known whether they share same universal bounds as the Bianchi-I model and if one 
required extra conditions to obtain the boundedness of these physical quantities. Further, 
the important role played by inverse triad modifications in Bianchi-IX spacetime was so 
far unrecognized. It is important to emphasize that the bounds on p, 9 and are one of 
the main distinguishing features between LQC and GR. In GR, the unboundedness of these 
quantities results in break down of geodesic evolution, a characetristic feature of singularities. 
In this work, we have established for the first time, that such divergences do not occur in 
the loop quantization of Bianchi-II and Bianchi-IX spacetimes in LQC. It has been shown 
for isotropic models, that these bounds play an important role for generic resolution of 
singularities in isotropic spacetimes in LQC |22|. Novel results obtained in this work, take 
us a step closer to prove generic resolution of singularities in Bianchi-II and Bianchi-IX 
spacetimes. 

Given the remarkable coherence of results for the isotropic and Bianchi-I models, and a 
promising strategy to quantize Bianchi-II and Bianchi-IX models which is consistent with the 
loop quantization of the Bianchi-I model, one may have expected that answers to the above 
questions would be on the lines of the results obtained in the previous models. However, we 
show that this expectation turns out to be not true and the answers leads to some surprises. 
We demonstrate that in the Bianchi-II model, the energy density has a global maxima. 
The expansion and shear scalars are bounded only if one imposes energy conditions on the 
matter content. In particular, it is important that the energy density has a lower bound 
(a simple requirement) which yields an upper bound on the shear scalar. The bounds are 
not universal and their values depend on the imposed energy conditions. Recall that for the 
isotropic and Bianchi-I models, the bounds on p, 6 and cx^ turned out to be universal and 
did not depend on energy conditions. It is to be noted that the in the effective Hamiltonian, 
matter Hamiltonian is treated as if Fock quantized. If matter is polymer quantized, in 
the way geometry is, results on dependence on energy conditions can a priori change. For 
the Bianchi-IX model, the inverse triad corrections are critical to obtain a bound on the 
energy density and the expansion and shear scalars when the approach to singularity in the 
classical theory is isotropic. In the Bianchi-IX model, energy conditions are shown to play 
little role. Inverse triad corrections are also important in other types of singularities, such 
as the cigar and pancake singularities, to achieve a finite value of the above quantities. In 
the effective description, singularity is resolved in the sense that the energy density, and the 
expansion and shear scalars are finite. However, there also exist following two mathematical 
possibilities. If the effective dynamics allows solutions where two of the triads tend to zero 
and the third diverges, then p^, 9 and o"^ diverge in Bianchi-IX spacetime. Also, if the 
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physical solutions exist such that one of the triads vanish and the other two diverge, then 9 
and (T^ diverge. We emphasize, that due to the underlying assumptions in the derivation of 
effective Hamiltonian constraint, it is not clear whether the effective dynamics description 
is valid when above two cases arise. 

Though the bounds on energy density, and the expansion and shear scalars for Bianchi-II 
and Bianchi-IX models turn out to be different from the Bianchi-I case, it is important 
to emphasize that these bounds are not optimal, as shown by a separate analysis (33| . 
Nevertheless, results on Bianchi-II and Bianchi-IX models stand in sharp contrast to those 
obtained for the isotropic and Bianchi-I models, where for the first time the role of energy 
conditions and the inverse volume modifications become important. Since the only change 
in the quantization strategy of these models, in comparison to the earlier ones, is in the way 
field strength is expressed in terms of connection, the cause of above differences lies in the 
usage of the connection operator approach. This quantization strategy is consistent with 



the holonomy based quantization of the Bianchi-I model [13| , but important differences can 
arise in presence of spatial curvature, as is evident from the analysis of Bianchi-IX model. 
The effective Hamiltonian of the loop quantization of Bianchi-IX spacetime does not lead to 



the one for the holonomy based loop quantization of = 1 model in the isotropic limit [10 
Rather, its isotropic limit corresponds to the connection operator quantization of isotropic 
spatially closed model [25|, where the expansion scalar does not have a maxima. We also 
find that connection operator approach for k = —1 model suffers from the same problem 
as shown in Appendix B. The unboundedness of 9 in the connection operator approach for 
the isotropic = ±1 models, does not directly affect its viability in the anisotropic models. 
However, these results show that this approach leads to a qualitatively different behavior, 
in comparison to the holonomy based quantization, in the presence of spatial curvature. 

We conclude with a discussion of some open questions resulting from our analysis. In 
this work, we have used the effective Hamiltonian framework for the entire range of vol- 
ume. In previous studies, extensive numerical simulations have confirmed the validity of 
effective spacetime description in isotropic models for universes where bounce occurs at vol- 



umes greater than Planck volume [2H6|, llO| and also Bianchi-I models (albeit for a different 
quantization) jss], 36|. It will be reasonable to expect the validity of bounds for the scales 
of interest, since the maximum values of physical quantities occur at volumes greater than 
the Planck volume. Nevertheless, this is an important issue which deserves a careful exami- 
nation in future. In particular to understand the behavior of energy density, expansion and 
shear scalars when approach to classical singularity is non-isotropic in the Bianchi-IX model. 
For this one has to derive effective Hamiltonian dynamics using embedding approach in the 
presence of anisotropics and examine its validity at small volumes. The pertinent question is 
whether such an analysis yields further corrections to the effective Hamiltonian constraint, 
and if so, how are the bounds affected in Bianchi-II and Bianchi-IX models. It will be inter- 
esting to see if such modifications bring the maximum values of physical quantities closer to 
the isotropic and Bianchi-I models. The second issue concerns with the energy conditions. In 
GR, these conditions play an important role in the proof of singularity theorems. Their role 
in LQC has so far been irrelevant. Analysis of modified Einstein's equations in Bianchi-II 
model reveal, that in LQC, expansion and shear scalars are unbounded if the energy density 
is unbounded below. This suggests that for cases where the latter is satisfied, singularities 



The inverse triad modifications lead to some interesting effects due to the choice of lapse in the Bianchi-IX 
model. These issues will be reported separately [sjj. 
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may not be resolved by quantum geometry effects in LQC Tliis feature confirms witli tlie 
general expectations that quantum gravity may not resolve all the singularities and weed out 
the unphysical situations with arbitrary negative energy, as considered above 37|. In future 



work, it will be important to investigate this issue in detail and classify the singularities 
corresponding to such solutions. Our analysis also shows the non-trivial role played by the 
inverse volume modifications in resolution of isotropic singularity in the case when spatially 
topology is compact. Perhaps, due to simplicity of the model, this was never revealed in the 
isotropic case.^^ However, for Bianchi-IX model it is the interplay of both the non-local field 
strength tensor and the inverse volume modifications which results in bounds on physical 
quantities. These results seem to point out that in general for spatially compact manifolds, 
physics of singularity resolution might be incomplete without inverse volume modifications. 
Finally, it is pertinent to revisit the connection operator approach, by demanding that the 
quantization of Bianchi-IX model be consistent with holonomy based quantization of /c = 1 
model. This is important as it will establish a consistency between the connection operator 
approach and holonomy based quantization in the presence of spatial curvature. It will be 
interesting to probe the way results obtained in this work are affected in such a quantization. 
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Appendix A: Expansion scalar in /c = 1 model 

In this appendix, we summarize the derivation of the expansion scalar using the effective 
Hamiltonian constraint of k = 1 model in LQC. Based on the techniques used for several 
models, to express field strength in terms of holonomies over a closed loop, the quantization 
of this model was performed in Ref. [lo|, which showed absence of singularity and a non- 
singular bounce when p = pmax~ 0.41ppi. More recently, an alternative quantization of this 



spacetime has been proposed [25| , which is based on the connection operator approach used 



for Bianchi-II and Bianchi-IX models [13|, |15|. Using the terminology from Ref. [25[, we will 
refer to these as "holonomy based" and "connection operator" quantizations respectively. 

The spatial manifold in = 1 model is S'^ with radius io = Vo^^, where the fiducial volume 
Vo = I6n^. In terms of the connection c and triad p variables, satisfying {c,p} = SnG'j/S, 
the classical Hamiltonian constraint is given by 



H 



cl 



3p2 



'2 



—7 

4 ' 



(Al) 



where e denotes the orientation of the triad. The triad p is related to the scale factor a{t) 
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Novel phenomenological implications have though been discussed, see for eg. 



Unlike the main body of this work, we will allow both the orientations of the triads in the appendices. 
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in the spacetime metric 



dr'^ 



1 — 



+ P {de^ + sin^ 



(A2) 



as p = ea'^i^. The equation of motion for the triad can be derived using the Hamilton's 
equations: 



P 



7 



et 



(A3) 



where the 'dot' refers to the derivative with respect to the proper time. From which one 
obtains the relation 

c = eio (^-fd + . (A4) 

Using these equations, it is straightforward to obtain the expansion scalar for the comov- 
ing observers: 6 = 3H = 3d/a, which turns out to be 



9 



7V bl 



(A5) 



In the classical theory, the approach to singularity is characterized by a divergence in the 
connection c, along with a divergence in 1/ \/\p\, causing the expansion scalar to become 
infinite. Let us now discuss the way this behavior changes in LQC for the holonomy based 
quantization performed in Ref. 10 and the connection based quantization (25| . 



Holonomy based quantization: For the loop quantization of = 1 model, proposed in Ref. 
lo| . the effective Hamiltonian constraint for lapse N = V is given by 



n 



eS 



3 



sm I /i I c — e 



X 



(A6) 



where x = sin^ (/^^t) " (1 + 7^)(Ai^|)^ /"^ = ^^/P and = iV^rr'jlp^. This includes the 
inverse triad modifications, which appear only in the matter Hamiltonian. The Hamilton's 
equation for triad yields 



2p 

p=-sin ^c-s^ 
and the expansion parameter turns out to be 



cos I /i I c — e 



9 = — sin ( /i ( c — e — cos ( jj,(c — e — 

This is a bounded function with a maximum value given by 

3 



9r, 



(A7) 



(A8) 



(A9) 



This is to facilitate a comparison between the isotropic hmits of the Bianchi model with the fc = 1 model. 
We thank E. Wilson-Ewing for a discussion on this issue in fc = 1 model. 
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which agrees with the maximum value in isotropic j22| and Bianchi-I model 31 1. It 



IS 



important to note that the above value is not affected by the inverse triad modifications, or 
the choice of energy conditions i.e. the specific form of matter Hamiltonian. In this sense, 
the maximum value of the expansion scalar is universal. 



Connection operator quantization: Based on the loop quantization of Bianchi-II and 
Bianchi-IX spacetimes, one can quantize the k = 1 model using an alternative method 
which avoids expressing field strength in terms of holonomies over a closed loop. Instead, 
the field strength is expressed in terms of the connection operator obtained from holonomies 
over open edges. This quantization results in the following effective Hamiltonian: 



n 



3p' 



eff 



87rG'72 



sin(/ic) 



/i 



^2 I + 



(AlO) 



The Hamilton's equation for the triad yields 



P 



7 



sin (2/ic) 



/i 



— 2e-^ cos (fic) 



(All) 



and the expansion parameter turns out to be 



27 



1 2 

— sin (2uc) — e— — -== cos (uc) 



(A12) 



If one assumes the validity of T-Les for all values of triad, the expansion scalar turns out to be 
unbounded: as p — )■ 0, ^ — )■ — oo.^^ This behavior is in agreement with that of the expansion 
scalar of the loop quantization of Bianchi-IX model without the inverse triad corrections. 

In conclusion, we find that there is an important qualitative difference in the behavior of 
expansion scalar in the holonomy based quantization (eq. flASp ) and the connection operator 
quantization (eq. flA12p ) for the k = 1 model. The expansion scalar in the holonomy based 
quantization turns out to be bounded by a universal value irrespective of the choice of 
matter, where as for the connection operator quantization it is divergent. Since the latter 
is tied to the existence of singularities, this result shows that with in the applicability of 
effective Hamiltonian constraint, connection operator quantization for k = 1 model has a 
drawback. 



Appendix B: Expansion scalar for k = —1 model in LQC 



Loop quantization of A; = — 1 model has been proposed in Refs. [ll|, Il2j. In this model, 
the Ashtekar-Barbero connection contains off-diagonal terms which poses technical problems 
in construction of the quantum theory using holonomies as elementary variables. In order to 



Effective dynamics for tliis model in the absence of tlie inverse triad corrections has been studied in Ref. 



25| , however problems with expansion scalar as found here were not discussed. 

It is possible that the effective Hamiltonian description may attain additional corrections as p — !■ which 
may regularize the above unbounded behavior. At the present stage, this is an open issue. 
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overcome these issues, loop quantization has been performed by considering holonomies of 
the extrinsic curvature. The quantum Hamiltonian constraint turns out to be non-singular 



with features similar to the quantization of other isotropic models Il2| . The resulting 



physics has been analyzed using the effective Hamiltonian constraint which reveals resolu- 
tion of strong singularities [l8[. A drawback of this quantization is that field strength is 
not expressed in terms of the holonomies of the connection^^, which motivates to consider 
alternate quantizations of k = —1 model in LQC A possible avenue would be to quantize it 
on the lines of Bianchi-II and Bianchi-IX models using connection operator approach which 
would bypass the above problem 39|]. Here we compute the expansion scalars in both of 
these approaches using the effective Hamiltonian constraint. As in the k = 1 model, we will 
refer to them as "holonomy based" and "connection operator" quantizations respectively. 
The classical Hamiltonian constraint for the k = —1 model (in lapse = ^) is given by 



3p^ 



c + e- 



1 



(Bl) 



where lo (not to be confused with £o introduced earlier) refers to the cube root of the fiducial 
volume Vo of the fiducial cell V needed to introduce the symplectic structure on the non- 
compact manifold. The triad p is related to scale factor as 



2/2 

p = ea 

and its time variation is given by 

. 2/H ( ^ 

Substituting the expression for p in terms of scale factors into the eq. IB3I we obtain, 

sL .J 

c + — = 'jealo , 
using which one finds the expansion scalar as: 



(B2) 



(B3) 



(B4) 



9 = 3- 



(B5) 



From the above expression, we find that the expansion scalar in the classical k = —1 model 
diverges as the singularity is approached. 

Holonomy based quantization: The effective Hamiltonian constraint for the loop quanti- 
zation of = — 1 model as performed in Refs. ll|, |l2| is given by 



H 



eff 



3 p^ 



sm I yU I c + £ 



+ X 



(B6) 



^•^ Nevertheless, it has been argued that in the symmetry reduced setting, due to gauge fixing it is possible 
to consider extrinsic curvature at the same footing as the connection 
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where x = "~7^(/^^^)^ ^^^d Ji is the same as in = 1 model. Since the spatial manifold is 
noncompact, there are no inverse triad corrections in T-LcK- A straightforward calculation, as 
performed for the k = 1 model shows that the expansion scalar turns out to be 



n 3 . / / lo\\ / / lo 

d = —^m \ u \ c + e— \ I cos \ u I c + e— 



(B7) 



where we have used the Hamilton's equation for the triad: 



2p . 



P = — sm ( /i ( c + e— 



cos ( /i ( c + e— 



(B8) 



The expansion scalar for the holonomy based quantization of = — 1 model thus turns out 
to be a bounded function with the same maximum value as the expansion scalar for A; = 1 
model feq. (lA9|) l 

Connection operator quantization: The effective Hamiltonian constraint in this case can 
be obtained by replacing c with sin(/ic)//i in the classical Hamiltonian constraint (IBip : 



eflf 



sin(/xc) lo_ 



The resulting Hamilton's equation for the triad becomes 



(B9) 



P 



7 



sin(2/ic) lo , ■ 
z 1- COS (/ic) 



(BIO) 



which leads to the following expression for the expansion parameter, 



3 
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A 



sin (2/ic) 



lo 2 



cos (/ic) 



(Bll) 



Unlike the behavior of the expansion scalar in the holonomy based quantization, Q does 
not have a maxima. The connection operator approach encounters the same limitation as 
the = 1 model regarding the behavior of the expansion scalar. Thus, there are important 
qualitative differences in the resulting physics from holonomy based and connection operator 
quantizations in LQC 
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